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Abstract 



The question addressed in this paper is: What originates pure spin su- 
perradiance in a polarized spin system placed inside a resonator? The term 
"pure" means that no initial coherence is imposed on spins, and its appearance 
manifests a purely self-organized collective effect. The consideration is based 
on a microscopic model with dipole spin interactions. An accurate solution 
of evolution equations is given. The results show that the resonator Nyquist 
noise does not play, contrary to the common belief, any role in starting spin 
superradiance, but the emergence of the latter is initiated by local spin fluc- 
tuations. The decisive role of nonsecular dipole interactions is stressed. 
A polarized spin system prepared in a nonequilibrium state returns to equilibrium 
through the spin-spin and spin-lattice relaxation mechanisms. The spin relaxation can 
be drastically accelerated if the nonequilibrium magnetic system is placed inside a coil of 
a resonance electric circuit with the natural frequency tuned to the precession frequency 
of spin magnetic moments [1]. The strong shortening of the relaxation time is caused by 
the coherence between individual rotating spins, which develops as a result of the interac- 
tion between the rotating magnetization and the resonator feedback field. This coherent 
phenomenon is analogous to the Dicke superradiance [2] occurring in atomic and molecu- 
lar systems, and so Bloembergen and Pound [1] also called this fast collective damping in 
spin systems the radiation damping. Priedberg and Hartmann [3] noted that, in fact, the 
whole process in spin systems involves no radiation at all, but merely nonradiative transfer 
of energy from the sample to the coil and back, the energy being dissipated in the circuit 
ohmically. And what one measures in experiments with spin systems is not the intensity 
of radiation, but the power of electric current. Nevertheless, the term superradiance has 
been accepted for the transient coherent phenomenon in spin systems, when it develops as 
a self-organized process, similarly to the Dicke superradiance. The similarity between the 
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latter and the spin superradiance is not the sole excuse for the accepted term. Another rea- 
son is to distinguish the spin superradiance from nuclear induction, free or collective, that 
is also a coherent phenomenon although not self-organized, but forced by thrusting onto a 
sample an initial coherence. One more justification for using the term spin superradiance 
is that this effect is always accompanied by coherent magneto-dipole radiation, though the 
corresponding radiation intensity is very weak as compared to the easily measured current 
power [4]. 

The term pure spin superradiance is used in order to stress that this is a purely self- 
organized process, when coherence develops from an absolutely incoherent state. This is to 
be distinguished from triggered spin superradiance during which collective effects also play an 
important role, but the process starts from a coherent initial state, so that the imposed initial 
coherence triggers the development of a correlated state, in the same way as triggered optical 
superradiance [5] happens in atomic and molecular systems. When spin superradiance is 
caused by nuclear spins, it can be called the nuclear spin superradiance. A system of ion 
spins in a resonator cavity can, in principle, be also a source of spin superradiance. 

The nuclear spin superradiance has been recently observed in a series of experiments [6- 
11] with different substances: from Al nuclear spins in ruby ( AI2O3 ) and from proton spins 
in propanediol ( C^H^02 ), butanol ( CiH^OH ), and ammonia ( NH^, ). The interpretation 
of the pure spin superradiance in these experiments was based on the following picture. A 
system of polarized spins is placed in a constant magnetic field directed opposite to the 
sample polarization. This means that the spin system is prepared in an inverted state. The 
sample is put inside the coil of a passive electric circuit whose natural frequency is tuned to 
the Zeeman frequency of spins. Fluctuating magnetic field formed by the thermal Nyquist 
noise of the resonance circuit starts moving spins from their position of unstable equilibrium. 
The motion of spins induces electric current in the circuit, which creates a stronger magnetic 
field acting back on spins. Under the action of the feedback field, spins move faster increasing 
even more the resonator feedback field, and so on. This avalanche-type process results in a 
fast spin relaxation. Such is the commonly accepted picture of pure spin superradiance. 
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However in this, generally correct, picture there is one suspicious point, namely, that the 
beginning of the process is originated by the thermal Nyquist noise of resonator. Really, 
if one attentively reads the classical paper by Bloembergen and Pound [1], then one finds 
there the estimate for the thermal damping, due to the thermal noise in resonator, showing 
that this damping is so negligibly small for macroscopic systems that it can never produce 
the initial thermal relaxation. 

Thus we confront the alternative: either the common belief that this is the resonator 
thermal noise which initiates the pure spin superradiance is a delusion or Bloembergen and 
Pound are wrong. To resolve this paradox and to answer the question "what actually is the 
origin of pure spin superradiance" is the aim of the present paper. 

The solution of the formulated problem meets the following difficulty. As follows from 
the analysis of Bloembergen and Pound [1], the homogeneous approach provided by the 
Bloch equation is not sufficient for correctly describing the process, but inhomogeneous lo- 
cal fields, that produce a microscopic relaxation mechanism, are essential. The phenomeno- 
logical Bloch equation, even being solved in a reasonably accurate approximation [12], can 
describe only the triggered spin superradiance, when an initial coherence is imposed by 
assuming nonzero initial conditions for transverse magnetization. To take into account in- 
homogeneous local fields providing a microscopic relaxation mechanism means the necessity 
of dealing with a microscopic model. Writing the equations of motion for spin components 
we get a 3A^ -dimensional system of nonlinear differential equations for N spins, plus the 
Kirchhoff equation for an electric circuit. If one invokes any approach, to solve this system 
of equations, based on the uniform mean-field approximation, then one immediately returns 
to a homogeneous picture equivalent to the Bloch equation, thus loosing an information on 
local fields. When the number of spins, N , is not too large, say iV ~ 10 — 10^ , then it is 
possible to resort to numerical calculations. Such a computer simulation has been done [4] 
(for mathematical details see [13]) and has been shown that, really, the pure spin superra- 
diance can exist without the resonator Nyquist noise. However, such computer simulations 
have the following deficiencies: (i) they are time consuming; (ii) they are able to give only 
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a qualitative description, since the number of spins involved is incomparably smaller than 
what one has in real samples with N ~ 10^^ ; (iii) they do not give analytical formulas 
that would be convenient to study with respect to the variation of all, sometimes numerous, 
parameters characterizing the system. Therefore, computer simulations can give a feel of 
what is happening, but cannot provide decisive answers. Below an analytical solution of 
microscopic equations is presented. 

A system of nuclear spins is described [14] by the Hamiltonian 

N N 

H^-^H,,-t,^BSi (1) 

iytj i=l 

with dipole spin interactions 



Hi^ — SiS-j 3 ( Sif^ii ] ( Sjf^i 



(2) 



where // is a nuclear magneton and 
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The total magnetic field 



B= Ho e, +H (3) 



consists of a constant external field Hq and an alternating field H of a resonator coil. 
The latter has n turns of a cross-section area A^. over a length / . The resonance electric 
circuit includes a resistance R , inductance L , and capacity C . The alternating resonator 
field 

is formed by an electric current satisfying the Kirchhoff equation 
in which Ef is an electromotive force and 



is a magnetic flux through the coil; rj = V/Vc being a flUing factor; Vc = lAc , the coil 
volume; p = N/V , density of spins. 



Define the resonator natural frequency lu = 1/y/LC , ringing width 73 = R/2L , and 
dimensionless fields 

h.f, /-^. (6) 

Introduce the parameter 

ao = nri— (7) 

characterizing the strength of coupling between the spin system and the resonator. Let us 
also use the notation 

1 TV 1 Af 

^^i^E(^r), «-^E(^n (8) 

i=i i=i 

for the mean spin components, where (. . .) implies statistical averaging. Then the Kirchhoff 
equation (5) takes the form 

dh /■* 

— + 273/1 + / h(r)dr = -2ao (u* + u)+ 73/. (9) 
at Jo 

To derive the evolution equation for the variables (8), we proceed as follows. Write the 
Heisenberg equations for the corresponding spin components with the standard notation 
cuq = iiHo/h for the Zeeman frequency. Decouple the double spin correlators in the manner 
described by ter Haar [15], in order to preserve the terms containing the homogeneous 
spin-spin relaxation 72 = , which can be done by using second-order perturbation 
theory. For generality, we may also include the term describing the spin-lattice relaxation 
7i = Tf" ^ . These steps are known and clear. The most difficult problem is how to treat the 
local spin fields 

1 3 
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caused by the dipole interactions 



^2 



r 



3^2 

bij = si"^ (-^2(/7jj) , (11) 



3„2 

= - jy sin^(2i?jj) exp {-i(pij) , 

where i^jj and ipij are the spherical angles of riij . Note that in a uniform approximation 
the local fluctuating fields (10) are zero because of the properties of the dipole interactions 
(11). To get a closed set of equations, at the same time retaining the information on the 
presence of fluctuating fields (10), we may replace the latter by stochastic fields, 5i — > 
^0, ^ ^ 1 the first of which, in compliance with (10), is real and the second is complex. 
The distribution of these random fields is such that the averaging over it, which we shall 
denote by ((...)), gives 

{{v^)) = m) = 0, {{^D) = \m')) = 7*^ (12) 

where the dispersion 7* , in accordance with (10), is of the order of 72 . In this way, for 
the spin components (8) we obtain the system of stochastic equations 

f = 2 ^^'^ + ^) - 2 ^^'^ + V^*) ^ - - 0, (13) 
= -272|m|^ - ^(73/1 + (p)su* + i(73/i + (p*)su. 

The structure of (13) is transparent: ^^h + cp is the total effective field acting on spins; h is 
the resonator field defined by (9); ipo and ip model random local fields with a distribution 
whose particular form is not important since all we need is the property (12). If ipo and 
(f were absent, then (13) would be reduced to the Bloch equation. 

To consider the case of pure spin superradiance, the initial conditions for the system of 
equations (9) and (13) are to be taken as 



h{0) = 0, u{0) = 0, s(0) = zq. (14) 

The electromotive force Ef — Eo cos u;t in (5) corresponds to the resonance mode of the 
thermal Nyquist noise. The driving force in (6) is 



cfiEi 



f^focosut; fo = ——^. (15) 

The system of equations (9) and (13) can be solved by a method [16] combining the 
guiding-center approach [17] and the method of averaging [18]. The idea is straightforward: 
First, we classify the variables onto fast and slow. To this end, we take into account the 
usual inequalities 71 <S a;o, 72 <S <^o, 73 <S c*; , and consider the quasiresonance case, when 
|A| <^ a;o , where A = uj — ujq is detuning. Thence we notice right away that the variables 
h and u can be treated as fast, while s and |Mp as slow. Keeping the latter as fixed 
parameters 

s — z, \u\ — V, (16) 

we get for the fast variables a system of linear equations, which, therefore, is not too difficult 
to solve. The found solutions for fast variables are to be substituted into the equations for 
slow variables, and the right-hand sides of these equations are to be averaged over the period 
27: /(Vq of fast oscillations and also over the distribution of stochastic fields. This procedure 
results in the equations 

^ = 972W - 7i(^ - C) - 7/^, 
dw 2 

— = -272W - 2g'^2wz + 27/^ (17) 

for the slow variables, where 

2 272 ^73^ 77(72-73)^ 

W = V 9 = OiQ 



,12 J (72 -73)^ + A2' 



and the attenuation 



in which 



(a-27r/?)(A'-72) + 272A(/? + 27ra) L (18) 



_ ,'73\ 7r(72 - 73)^ a (lA 7^(72 -73) A 



.a;o/ (72 - 73)^ + A^ ' \ujqJ (72 - 73)^ + A^ ' 

is due to the action of the driving field (15). 

The amphtude of the electromotive force related to the thermal Nyquist noise [19], at 
temperature T satisfying the inequahty hw <^ ksT , is given by Eq — ^^RkBT/n . 
Whence for the amplitude of the driving field (15) we have 

^° - ^^^^ 

Let us accept the values of parameters characteristic of experiments [7-11] with proton 
spins: a;o ~ a; ~ 10^ sec~^, 71 ~ 10"^sec~^, 72 ~ lO^sec"^ 73 ~ lO^sec"^, T ~ O.IK , 
and N ~ lO^^ ^hen /o ~ 10"^° and the thermal attenuation (18) is 7/ ~ lO-^^sec"^ . 
Such an insignificant quantity, of course, plays no role, as compared to all other damping 
parameters, and has to be neglected in (17). 

This result shows, in agreement with Bloembergen and Pound [1] , that the Nyquist noise 
of resonator can never produce the initial thermal relaxation, thus, cannot originate the pure 
spin superradiance. 

Ommitting in (17) the negligibly small 7/ and taking into account that 71 <^ 72 , we 
come to 

dz dw , . 

_ = ^72U), — = -272U)(1 + gz). (20) 

According to (14), the initial conditions are z{0) — zq and v{0) = . Equations in (20) 
are exactly integrable yielding 

z = tanh 



9l2 V To / g 

2 



>L rfC^UH:|.; (21) 

V^72y V To y uji 



here 70 is the radiation width given by 



7o = - 2(t/72)^£*^o, 



(22) 



where 



To = 72(1 + 5'2;o), 



oJqJ 



the radiation time tq = 70 ^ ; and the delay time is 



to = -ln 



7o - To 



(23) 



7o + To 

The criterion for the occurrence of spin superradiance is the validity of the inequalities 



< < CO, Tq<T2 



(24) 



Invoking (22) and (23) and bearing in mind that -C 1 , we find that (24) is equivalent to 



Zq < Zp = , 
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£* > 0. 



(25) 



As far as \zq\ < 1/2 , the first of the inequalities (25) requires that g > 4: . In this way, the 
pure spin superradiance occurs when the initial spin polarization zq is negative, with an 
absolute value surpassing the threshold \zp\ = 2/g , when the coupling of the spin system 
with a resonator is sufficiently strong, g > i , and if there exist local random fields with a 
nonzero dispersion 7* > . 

To emphasize the crucial importance of the local fields, let us notice that if one puts 
£* ^ , then 7o —>■ \To\ and |to| —>■ 00 . That is, without these fields the pure spin 
superradiance is impossible. To make the essential dependence of the delay time on 
apparent, we may write (23) for the case of strong coupling, when g\zo\ ^ 1 , then 



to 



To 



'^g\zo 



In 



2zn 



Prom here it is evident that to ^ 00 as ^ . For the parameters typical of experiments 
[7-11] we have to ^ 10~^ — lO^^sec . So, this is the local random fields that are responsible 
for starting the process of self-organization leading to the pure spin superradiance. 
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One more question is worth answering: Which part of the local fields is more important 
for initiating the pure spin superradiance? Recall that the stochastic fields entering into the 
evolution equations (13) are related to two types of local fields defined in (10). As follows 
from (13), the term 5i in (10), corresponding to (fo , only shifts the rotation frequency, 
while the term (pi , corresponding to (p , starts moving the spin z -component even when 
the resonator feed back field h is yet absent. The term ipi in (10) is due to the dipole 
interactions bij and defined in (11). These interactions, in the theory of magnetic 
resonance [14], are called nonsecular interactions, as compared to the secular interaction 
Gij . The initial motion of spins, when u{0) = and h{0) = , is due solely to the action 
of nonsecular interactions. This conclusion is in agreement with computer simulations [4,20] 
for small and mesoscopic spin systems with ~ 10 — 10^ . Thus, we are in a position to 
give the final answer to the question posed in this paper: 

The pure spin superradiance in a nonequilibrium system of polarized nuclear spins can 
be originated only by local fields due to nonsecular dipole interactions. The thermal Nyquist 
noise of resonator plays no role in this process. 
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